Energy Spectrum and Quantum Hall Effect in Twisted Bilayer Graphene 
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We investigate the electronic structure and the quantum Hall effect in twisted bilayer graphenes 
with various rotation angles in the presence of magnetic field. Using a low-energy approximation 
which incorporates the rigorous interlayer interaction, we computed the energy spectrum and the 
quantized Hall conductivity in a wide range of magnetic field from the semi-classical regime to the 
fractal spectrum regime. In weak magnetic fields, the low-energy conduction band is quantized into 
electron-like Landau levels and hole-like Landau levels at energies below and above the van Hove 
singularity, respectively, and the Hall conductivity sharply drops from positive to negative when the 
Fermi energy goes through the transition point. In increasing magnetic field, the spectrum gradu- 
ally evolves into a fractal band structure called Hofstadter's butterfly, where the Hall conductivity 
exhibits a non-monotonic behavior as a function of Fermi energy. The typical electron density and 
magnetic field amplitude characterizing the spectrum monotonically decrease as the rotation angle 
is reduced, indicating that the rich electronic structure may be observed in a moderate condition. 
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I. INTRODUCTION 

The electronic structure of bilayer graphene is highly 
sensitive to the stacking geometry between the two lay- 
ers. The interlayer interaction in bilayer graphene with 
regular AB stacking^— modifies the linear dispersion of 
monolayer graphene into the quadratic dispersion, where 
an electron behaves as a massive particle* 4 - On the other 
hand, the recent epitaxial growth technique^ realized 
twisted bilayer graphene (TBG) in which two layers are 
stacked with a random rotation angled— The unit cell 
area of TBG can be more than a thousand times as large 
as that of monolayer graphene, due to slightly misori- 
cntcd lattice vectors of two layers. Such an atomic con- 
figuration was observed as Moire pattern in the scanning 
tunneling microscopy^— TBG was also fabricated in dif- 
ferent methods such as folding of mechanically exfoliated 
graphenesjM segregation of graphene on Ni filing and 
unzipping of carbon nanotubeJ^ 

The electronic structure of TBG shows a linear band 
dispersion near Dirac pointsii 1 ^— rather than the mas- 
sive dispersion of AB stacked bilayer, suggesting rela- 
tively weak interlayer interaction. In strong magnetic 
fields, however, it is predicted that the spectrum ex- 
hibits a fractal structure called Hofstadter's butterfly, in 
which a series of energy gaps appears in a self-similar 
fashion . 20 ' 21 The fractal band structure generally occurs 
in a periodic system when the magnetic flux per a unit 
cell is comparable to h/e, and this condition is real- 
ized in TBG in a reasonable magnetic-field range owing 
to the large unit cell. The fractal band structure and 
the quantum Hall effect were theoretically studied for 
TBG in the strong magnetic field regime using a con- 
tinuous interlayer coupling modeli^i Experimentally, the 
energy spectrum of the twisted graphene stacks in mag- 
netic field was probed in the transport measuremen t 22 ' 24 
and the magneto-optical absorption^ while the fractal 
band structure has not yet been observed. 

In this paper, we investigate the electronic spectrum 
and the quantum Hall effect in TBG with various rotation 



angles and magnetic fields. We calculate the spectrum by 
including a limited number of bases which are significant 
in the low-energy spectrum, while rigorously taking ac- 
count of transfer integrals between lattice points on the 
different layers. Using this method, we describe the spec- 
tral evolution in a wide range of magnetic field, from the 
semiclassical Landau levels in the weak field regime to 
the fractal band structure in the strong field regime. 

In weak magnetic fields, we find that the low-energy 
conduction band is quantized into electron-like Landau 
levels and hole-like Landau levels at energies below and 
above the van Hove singularity, respectively, in accor- 
dance with the topological change of the Fermi sur- 
face from electron-type to hole-type at the band saddle 
point. As a consequence, the quantized Hall conductiv- 
ity abruptly jumps from positive to negative when the 
Fermi energy goes through the transition point. In in- 
creasing magnetic field, the electron and hole Landau 
levels begin to be mixed and gradually evolve into the 
fractal band structure. We calculate the quantized Hall 
conductivity for each single gap, and demonstrate that it 
changes non-monotonically as a function of Fermi energy 
and magnetic fieldj^^ 



II. THEORETICAL METHODS 

A. Atomic structure 

TBG is characterized by the relative rotation angle 9 
and the relative translation vector between two graphene 
layers. When the lattice structures of the two layers are 
commensurate, we can define the primitive lattice vec- 
tors Li and L2 as the least common multiples of the 
unit vectors on the two layers. Li is written by integers 
m, n, to', n' as 2 ^ 
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where a^ and a 2 ' ' are the lattice vectors of the layer 
I = 1, 2 defined in Fig. [Ha). L2 is obtained by rotating 
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Li by 60°. By appropriate choice of lattice vectors af\ 
the indices (m', n') can be made equal to (n, m), and thus 
TBG is specified by a single pair of integers (m, n). The 
rotation angle 9 is related to (m, n) by 



1 TO 2 

cose/ = -- 



4mn 



2 m 2 + rt 2 + n 
and the lattice constant of TBG by 



L = |Li| = a^™ 2 + n 2 + 



mn, 
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where a = |ai| = |aa| ps 0.234 nm is the lattice constant 
of monolayer graphene. The area of TBG unit cell is 
given by S = |Li x L 2 | = (V3/2)L 2 . 

Figure Ufa) shows the atomic structure of TBG with 
(m, n) = (1, 2) and 9 = 21.8°. Throughout the paper, we 
set the coordinates (x, y) on graphene plane so that y axis 
is parallel to L 2 , and z to the direction perpendicular to 
the plane. We ignore the relative translation between two 
layers, which makes a minor difference in the electronic 
structure when the unit cell is large enough. 

Figure QJb) shows the extended Brillouin zone of TBG 
with 9 — 21.8°. The two large hexagons represent the 
first Brillouin zones of layer 1 and 2, respectively. K^ 1 ' 
and K'^ denote the two inequivalent corners of layer I, 
which are Dirac points in the single-layer band structure. 
The four Dirac points of K' (1 \ , and K'^ are 
folded back to two Dirac points, K and K', in the reduced 
Brillouin zone.— 

Figure [5] shows the atomic structures of four different 
TBGs to be considered in following sections. They are 
specified by (m,n) = (3, 4), (8, 9), (12, 13), and (22,23), 
and the rotation angles 9 = 9.43°, 3.89°, 2.65°, and 1.47°, 
respectively. As the angle 9 decreases, the size of the unit 
cell enlarges and the Moire pattern becomes evident. 
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B. Tight-binding model 

In a tight-binding model in terms of p z atomic orbitals, 
the Hamiltonian of TBG at zero magnetic field is written 

as 



ttB=0 

-"tbg 



£*(Ri,H,-)l*<><*il+hx.. (4) 

{id) 



d ■ e. 



where R^ and |\E^) represent the lattice point and the 
atomic state at site i, respectively, and t(Ri,Rj) is the 
transfer integral between the sites i and j. We adopt an 
approximation;^— 

/•:R,.R, : l ; , ; , T 1— ( — ; — ) +V rll ,T 

VppTT = Vppu eX P 

Vppa = V° pa cxp (- ° ) , (5) 
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FIG. 1: (Color online) (a) Atomic structure of TBG with 
rotation angle 6 = 21.8°. Dashed (red) and solid (green) 
lines represent the lattices of layer 1 and 2, respectively, (b) 
Brillouin zone of TBG with 9 = 21.8°. Dashed (red) and solid 
(green) large hexagons correspond to the first Brillouin zone 
of layer 1 and 2, respectively, and thick small-hexagon to the 
reduced Brillouin zone of TBG. Open and filled circles are 
two inequivalent valleys K and K' of TBG. 



where d = R; — Rj, and e z is the unit vector parallel 
to z axis. Vpp W is the transfer integral between the the 
nearest-neighbor atoms of monolayer graphene which are 
located at distance a = a/V3 w 0.142 nm, and Vp pa is 
the interlayer transfer integral between vertically located 
atoms at the interlayer distance do ps 0.335 nm. Here 
we take « -2.7 eV, V^ pa « 0.48 eV, to fit the low- 
energy dispersion of bulk graphite. 8 is the decay length 
of the transfer integral, and is chosen as 0.184a so that the 
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next nearest intralayer coupling becomes 0.1V 
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(m,n) = (3,4), L = 1.50nm 



(b) = 3.89° 




(m,n) = (8,9), L = 3.62 nm 



(c) = 2.65° 




(m,n) = (12,13), £. = 5.33 nm 



(d) 0= 1.47° 




(m,n) = (22,23), L= 9.59 nm 



FIG. 2: (Color online) Atomic structure of TBG with rotation angles of (a) 9.43°, (b) 3.89°, (c) 2.65°, and (d) 1.47°. Dashed 
(red) and solid (green) lines represent lattices of layer 1 and 2, respectively, (m, n) is the index characterizing the primitive 
lattice vector of TBG, and L is the length of the lattice vector. 



The transfer integral for d > 4ao is exponentially small 
and can be safely neglected. The band velocity of the 
Dirac cone in monolayer graphene is given by 



V <-> v rm- 



(6) 



We plot the band structures of four TBGs with the 
different rotation angles in Fig. [31 The band structure 
near the Dirac point exhibits a linear dispersion analo- 
gous to monolayer's,^^— while the interlayer interac- 
tion strongly influences the band structure in the higher 
energy region, giving the van Hove singularity to M 
poin t 10 ' 19 ' 33 ' 34 and a hole-like pocket to V point. As 
the rotation angle decreases, the energy span between 
the saddle point and Dirac point is shortened, and the 
Fermi velocity decreases accordingly] 8 ' 14 i 17 In the small- 
est rotation angle = 1.47°, in particular, the lowest 
energy band is highly distorted and nearly flat disper- 
sion appears near zero energy . 32 ' 35 A different choice of 
parameters slightly modifies the band structure of TBG 
but does not make qualitative differences. 



C. Electronic structures in magnetic fields 

We consider TBG in a uniform magnetic field B — 
(0, 0, B) perpendicular to the layer. For simplicity, we 
neglect spin Zeeman splitting throughout the paper. The 
system is characterized by the number of magnetic flux 
per a unit cell, $ = BS, measured in units of the flux 
quantum $o = h/e. In magnetic field, the Hamiltonian is 
no longer translationally-symmetric because of the spa- 
tial dependence of the vector potential. When $/<I>o is a 
rational number p/q (p and q are coprime integers), how- 
ever, we can introduce a magnetic unit cell with lattice 



vectors Li = qJj± and L2 = L2, and construct the eigen- 
states so as to satisfy the magnetic Bloch conditionj 36 i 37 
By choosing the vector potential as A = (0, Bx, 0) and 
taking y axis parallel to L2 , the magnetic Bloch condition 
for TBG is written as 



# k (r + Li) = e ik-L le - i (e/ft)(A-Bxr).L 1 ^ k ^ 



tf k (r + L 2 ) 



ik-L 5 



*k(r), 



(7) 



where k is the Bloch wavenumber defined in the magnetic 
Brillouin zone spanned by reciprocal vectors of Li and 
1*2- Since the magnetic unit cell is q times as large as the 
unit cell in the absence of magnetic field, the magnetic 
Brillouin zone is g-fold of the original, and each energy 
band at zero magnetic field splits into q subbands^ . 

The tight-binding Hamiltonian under a magnetic field 
is obtained by adding a phase factor to the transfer inte- 
gral in Eq. This is written as 

Htbg = ~J2 *( R *> R i) e#w l*i><*jl + h - c -> 



R; 



A(r) ■ dr. 



(8) 



Ri 



It is, however, not practical to calculate the energy spec- 
trum of TBG by diagonalizing this Hamiltonian, since 
the number of atoms in a magnetic unit cell is huge in 
feasible magnetic fields. Instead, we construct the basis 
from the effective mass wavefunctions for Landau levels of 
monolayer graphene, which approximate the eigenstates 
in the absence of the interlayer coupling. We then trun- 
cate the bases far from the Dirac point, and compose 
the Hamiltonian matrix by writing -Htbg m terms of the 
reduced basis. 

In monolayer graphene under magnetic field, the eigen- 
states are labeled by (v, n, k y ) with the valley index 
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(a) 0=9.43° (b) = 3.89° (c) = 2.65° (d) 0=1.47° 
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FIG. 3: Band structure of TBG with rotation angles of (a) 9.43°, (b) 3.89°, (c) 2.65°, and (d) 1.47°, in the absence of magnetic 
field. Dirac point energy is set to zero. 



v = K,K', the Landau level index n = 0, ±1,..., and 
the wave vector k y along y direction;^— The eigenen- 
ergy depends only on n as 



hw B sgn(n) v / |n|, 



with hhJB = V2hv 2 eB. The effective wavefunctions are 
written as 



,39,40 



F Knky (r) = 



( sgn(n)(-i)0| n |_i^ H (a;) ^ 

<t>\n\,k y (x) 




K'nk 



,(r) 
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a_ e ik yy 
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o 



<t>\n\,k y {x) 

\ sgn(n)(-i)<j>\ n \-i ik (x) 



(.10) 



Here F = (Fj£ , F§ , F% , F§ ) is a four-component 
vector representing the envelope function of 
each site and valley. We defined 4> n ,k{x) = 

(2 n n\^l B y 1/2 e- z2 / 2 H n (z), with z = (x + kl 2 B )/l B 
and H n being the Hermite polynomial, Ib — yh/JeB), 
and 



sgn(n) 



1 (n = 0), 

1/V2 (n^O), 

(n = 0), 

n/\n\ (n^O). 



(11) 



The tight-binding wavefunction ^ on the layer / can 
be expressed in terms of the envelope function F as^ 



^a(Ra) 
$b(Rb) 



FZiR^ + e^j^^Ff (R, 



- we l " <!) e lK(!, - RB Ff(R B ) 



+^ K ' (!) - Rs Ff (K B ), (12) 



where 77 w is the angle of a,'' to x axis. We define 

' ^ -L UlhfXry 

as the tight-binding wavefunction on the layer Z generated 
from F vnky . 

We then combine the bases of different k y so as to 
(9) satisfy the magnetic Bloch condition, Eq. ([7]). We define 



(0 



E 



or exp 



iirpq- 



j(j + 1) 



(0 



exp [i(k -^.(Li-if,; 

^-(K^-^fo'+m), (13) 

Ivy 



where k is the Bloch wavenumber defined in the mag- 
netic Brillouin zone, m = 0, 1, • • ■£>— 1, and represent 
KW,K'W for u = K,K', respectively. It is straightfor- 
ward to show that this satisfies the condition of Eq. (J7J . 

An eigenstate of TBG is written as a linear combi- 
nation of single-layer eigenstates ^fl\, m belonging to 
the same k. We only include single-layer bases within 
-£ max < En < Emax, to discard the bases which do not 
much affect the low-energy spectrum. The eigenenergies 
are obtained by diagonalizing the Hamiltonian matrix 
within the reduced bases, 



H k [(l,v,n,m),(l',v',n',m')} = (^IJ^tbgI*^ 



(14) 



for each k in the magnetic Brillouin zone. The cut-off 
energy should be sufficiently larger than the interlayer- 
coupling energy, which is of the order of Vp pa . at most, 
and tends to decrease in small twisting angles. Here we 
take E max = 1.5 eV for 6 = 9.43° and 3.89°, and 1.0 eV 
for 2.65° and 1.47°. To avoid undesired effects caused 
by a discrete change in the number of bases in varying 
magnetic field, we adopt a soft cut-off which gradually 
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reduces the matrix elements associated to the single-layer 
bases beyond ±E max . 

We calculate the matrix elements [Eq. (fl"4"l) ] between 
different layers by evaluating the transfer integral for each 
pair of carbon atoms up to the cut-off distance d = Aciq. 
The matrix elements within the same layer can be re- 
placed with a diagonal matrix composed of the effective- 
mass eigenenergies in monolayer graphene, 

i/k [(I, v, n, m), (I, v', n', m')] (15) 

This treatment is valid in low energies, as long as the 
magnetic field is not too strong, or lg 3> a. 

When the Fermi energy Ep is inside a band gap of the 
spectrum, the Hall conductivity <j xy is evaluated by the 
formul a 43 ' 44 



"'mL' 



where np is the electron density per unit area below the 
gap. 

III. RESULTS AND DISCUSSION 

We show the energy spectrum (left) and quantized Hall 
conductivity (middle) against the magnetic field ampli- 
tude, for 9 = 9.43°, 3.89° in Fig. H and for 2.65°, 
1.47° in Fig. [5] In the right-most panel, we show 
the zero-field band structure in the same energy range. 
The energy spectrum of 9 — 9.43° [Fig. Ufa)] is almost 
equivalent to monolayer's Landau level, suggesting that 
two layers are nearly decoupled in this energy region. 
The sequence of the Hall conductivity, 4, 12, 20, • • • in 
units of — e 2 //t ) 24 ' 45 ' 46 is exactly twice as large as the 
monolayer's . 1 ' 2 ' 40 ' 41 Each Landau level is eight-fold de- 
generate due to the number of layers as well as the spin 
and valley degeneracies. 

In contrast, the energy spectrum of 9 = 3.89° [Figs. 
BJd) and (He)] exhibits a complicated structure which is 
clearly distinguished from monolayer graphene. In weak 
magnetic fields of < 0.1, the low-energy spectrum 

below 0.2 eV shows monolayer-like Landau levels and Hall 
conductivity of 4, 12, 20, • • • . In the higher energy re- 
gion above 0.2 eV, on the other hand, we observe hole- 
like Landau levels moving downward in energy, and the 
negative Hall conductivity of 0, —4, —8, —12, • • • . When 
the electron density increases from the charge neutral- 
ity point, the Hall conductivity rises in a sequence of 
4, 12, 20, • • • with a step of 8, then abruptly drops to a 
negative extremum, and increases with a step of 4 all the 
way to zero. 

Those spectral features in weak magnetic field per- 
fectly coincide with the zero-field band structure in Fig. 
0Jf). The electron-like Landau levels are regarded as the 
quantized orbits accommodated in electron pockets at 
K and K' points, while the hole-like Landau levels are 
those in a hole pocket at Y point. The transition from 



electron-like levels to hole- like levels corresponds to topo- 
logical change of the Fermi surface at the saddle point 
(M point), which is responsible for the van Hove singu- 
larity at 0.2 eV. The step of the Hall conductivity reflects 
the number of electron and hole pockets in the first Bril- 
louin zone, i.e., the degeneracy of an electron- like level is 
twice as large as that of a hole-like level, because there 
are inequivalent K and K' points whereas there is only 
one r point. Note that the lowest conduction band is 
composed of a pair of nearly degenerate energy bands 
as shown in Fig. E{f), and this doubling contributes to 
the degeneracy of two in addition to the spin degener- 
acy. These two bands are not completely degenerate but 
mirror-symmetric with respect to the lines of K — Y, 
K' — Y and K — K' , and thus give the identical Lan- 
dau level energies. The corresponding states of the two 
bands are related by Ci symmetry in the real-space lat- 
tice structure with respect to the sides of the unit cell in 

Fig.m 

The electron-like and hole-like Landau levels are al- 
ternatively explained by a nearly-free electron model, 
without mentioning the rigorous zero-field band struc- 
ture. In Fig. [SJ we illustrate semi-classical electron tra- 
jectories at several different Fermi energies for a "free" 
TBG with interlayer coupling neglected. In the limit of 
a small Fermi energy [Fig. [H^a)], electrons move along 
closed orbits around K and K', and those motions are 
quantized into monolayer-like Landau levels. Since each 
of K and K' points include two original K points from 
top and bottom monolayers, the Hall conductivity yields 
4, 12,20, i.e., double of monolayer's sequence. For 
large Fermi energies, the electron orbits around the K 
and K' valleys cross each other as shown in Fig. [Sic). A 
finite interlayer coupling interchanges the orbits at each 
crossing point, and generates a single hole-like trajectory 
moving around Y point in the opposite direction. The 
corresponding hole- like Landau levels are four-fold degen- 
erate due to spin and the Fermi circle doubling, and thus 
the Hall conductivity takes 0, —4, —8, —12, ■ ■ • . The mid- 
dle panel [Fig.[6jb)] is for the intermediate energy region 
between two regimes. There, the different semiclassical 
orbits are strongly mixed by the magnetic breakdown due 
to a small fc-space separation, resulting in broadening of 
Landau levels near the van Hove singularity in Fig.^d). 

The electron density to fill the lowest conduction band 
is given by 

* 9 ' (17) 



n 



S 



where g s is the spin degeneracy and 2 is the band dou- 
bling, no characterizes the order of the electron density 
required to reach the van Hove singularity and the hole- 
like Landau levels. We have no = 3.5, 1.6 and 0.5 in units 
of 10 13 cm- 2 for 9 = 3.89°, 2.65° and 1.47°, respectively. 
The folded band structure and the corresponding Lan- 
dau level spectrum of TBG are quite analogous to the 
whole 7r band in monolayer graphene^, while the elec- 
tron density to access the hole-like levels is of the order 
of 10 15 cm" 2 . 
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FIG. 4: (Color online) Energy spectrum and quantum Hall effect in TBG in magnetic field with rotation angles of 9.43° (above) 
and 3.89° (below). In each row, the left and middle panels display the energy spectrum and the quantized Hall conductivity 
as functions of magnetic field strength, respectively, and the right panel shows the band structure in the absence of magnetic 
field. In (b) and (e), the quantized values of Hall conductivity inside energy gaps are indicated by numbers as well as colors 
filling the gaps. The Hall conductivity of the gray area cannot be determined by the present calculation. 



When the magnetic held amplitude becomes as large as 
$/$o ~ lj the uncertainty in electron momentum is com- 
parable to the size of the Brillouin zone, and the semi- 
classical picture breaks down. There the energy spectrum 
exhibits a fractal band structure ! 20 ! 21 In Fig. [4j we ob- 
serve that the electron and hole Landau levels gradually 
evolve into the fractal structure as <£/$o approaches 0.5, 
and the Hall conductivity in the fractal regime behaves 
non-monotonically as a function of Fermi energy . 25 ' 26 

The energy spectra of TBG with 9 — 2.65° and 



1.47° [Fig. [5] exhibit similar fractal structures to that 
of 9 = 3.89°, as expected from the resemblance between 
the band structures in the absence of magnetic field. The 
held amplitude corresponding to the same $/<5>o is lower 
for smaller 9, so the actual magnetic field strength needed 
to observe a fractal band structure becomes more feasi- 
ble in smaller rotation angle. The fractal feature be- 
comes conspicuous roughly when $/$o ~ 0.3, which cor- 
responds to 100 T, 50 T, and 15Tin6> = 3.89°, 2.65°, and 
1.47°, respectively. 
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FIG. 5: (Color online) Plots similar to Fig. [4] for TBG with rotation angles of 2.65° (above) and 1.47° (below). 




FIG. 6: Fermi circle and electron trajectories of TBG in a 
nearly free electron picture, for three different Fermi energies 
(a) in the vicinity of Dirac points, (b) near van Hove singu- 
larity at the saddle point, and (c) hole-like band at F point. 



In the case of 9 = 1.47°, the spectrum is strongly com- 
pressed in the vicinity of Dirac points, in accordance with 
the band width reduction which occurs in small rotation 
angles . 19 ' 32 ' 35 Although the band structure near Dirac 
points is almost flat, T point still has a finite band ve- 
locity. As a consequence, the energy gaps between the 
hole-like Landau levels are much wider than those be- 
tween the electron-like levels. 
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IV. CONCLUSION 

We investigated the electronic structure and the quan- 
tum Hall effect in TBG with various rotation angles in 
the presence of magnetic field. We calculated the en- 
ergy spectrum and quantized Hall conductivity in a wide 
magnetic-field range, and described the evolution from 
the semi-classical Landau levels to the fractal band struc- 
ture. In weak magnetic field, the low-energy conduc- 
tion band is quantized into electron-like Landau levels 
and hole-like Landau levels in accordance with the struc- 
ture of the folded energy band. In increasing magnetic 
field, those semiclassical levels gradually evolve into Hof- 
stadter's butterfly, where the Hall conductivity exhibits 
a non-monotonic behavior as a function of Fermi energy. 
The typical electron density and magnetic field ampli- 
tude characterizing the spectrum monotonically decrease 
as the rotation angle is reduced, indicating that the rich 



electronic properties may be observed in a moderate con- 
dition for TBG with small angle less than 5°. 
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